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Numerical  Test  of  Kelvin's  Phase  Approximation  in  Calculating 
the  Magnetic  Field  of  an  Elemental  Ring 
of  Meridional  Currents  on  a Sphere 

1.  INTRODUCTION 

The  numerical  error  introduced  by  using  Kelvin's  stationary  phase  approxima- 
tion in  the  calculation  of  the  magnetic  field  produced  by  an  elemental  ring  of  meri- 
dional current  moments  is  examined  in  this  report.  Such  a calculation  arises,  for 
example,  in  predicting  the  phase  of  the  1.1  I.OKAN  signal  propagating  over  a spheri- 
cal earth. 

An  exact  expression  for  the  magnetic  field  is  derived  in  the  form  of  a complex 
integral  over  the  entire  ring.  This  expression  then  is  reduced  via  Kelvin's  approxi- 
mation to  only  two  contributions:  one  from  the  nearest,  the  other  from  the  farthest 
current  moments.  The  complex  integral  is  numerically  integrated  via  Chebyshev's 
quadratures  and  compared  with  the  Kelvin  results. 

2.  THEORETICAL  ANALYSIS 

Assume  a ring  of  width  A0(m)  at  the  intersection  of  the  sphere  r = a with  the 
cone  0=  constant  and  containing  a uniform  meridional  surface  current  density 
jg  (A/m).  The  directions  of  the  current  moments  are  indicated  in  f igure  1 by  arrows 

(Received  for  publication  7 April  1978) 

1.  Wlntner,  A.  (1945)  Remarks  on  the  method  of  stationary  phases,  J.  Math.  Phys. 
24:127-130. 
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tangent  to  the  sphere  along  the  ring  AQ'QHA  and  pointing  toward  the  pole  T.  Let 
the  observation  point  P lie  on  the  sphere  at  the  polar  angle  8 ( in  the  meridian  plane 
which  passes  through  the  Cartesian  x - axis.  Then 

P (x,  y,  z)  = P (a  sin  8 Q,  0,  a cos  8 ) . 

A point  Q on  the  moment  ring  has  coordinates 

Q(x,  y,  z)  = Q (a  sin  8 cos  0,  a sin  8 sin  0.  a cos  8 ) . 

Consequently,  the  magnetic  field  produced  at  P by  the  oscillating  current  moment 
at  Q is  given  by 


A i a sin  9 da  ...  ... 

■ 9 y exp  (lkH) 

477  R 


dH(P)  = — ££P(iklil-  <-ik  + 4>(0  x R)  . 


(1) 


omitting  the  spinor  exp  (-iwt).  Here  k = 2 n / \ = free  space  wavenumber; 

R = a \[2  (1  - cos  8 cos  8 - sin  8 sin  8 cos  0) 1 ^ 2 
= distance  along  QP; 

R = ^ [x  (sin  8 cos  0 - sin  8 Q)  + y sin  8 sin  0 + z (cos  8 - cos  8 o>] 

= unit  vector  along  QP; 

8 - x cos  8 cos  0 + y cos  8 sin  0 - z sin  8 
= unit  vector  along  the  current  at  Q; 

and  x,  y,  and  z are  Cartesian  unit  vectors. 

The  magnetic  fields  produced  by  the  moments  at  each  pair  of  conjugate  points 
Q (8,0)  and  Q ' (8  , -0)  cancel  in  the  radial  and  meridional  directions.  Thus,  the  total 
magnetic  field  at  P is  strictly  azimuthal  and  given  by 


2.  Stratton,  J.A.  (1941)  electromagnetic  Theory.  McGraw-Hill,  New  York  and 
London,  p 435. 
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Figure  1.  Meridional  Current  Moments 
on  a Spherical  Surface  W ltd  Cartesian 
Coordinate  Axes  at  Its  Center 
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2u 

I a sin  8 
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(-  ik  + -jy ) 


• a (cos  0 - sin  8 sin  8 Q - cos  8 cos  8 cos  0)  . 


(2) 


Note  that  the  phase  factor  kH  oscillates  rapidly  when  ka  is  a large  number,  but 
that  the  rest  of  the  integrand  is  smoothly  varying  in  0.  By  Kelvin's  principle  the 
integral  then  need  be  evaluated  only  in  the  vicinities  of  stationary  points  in  the  phase 
factor.  Furthermore,  the  smoothly  varying  part  of  the  integrand  may  be  treated 
as  if  it  were  constant  near  any  one  stationary  point  and  factored  out  of  the  integral. 
The  phase  factor  itself  then  may  be  expanded  to  second  order  about  this  particular 
stationary  point  and  integrated.  This  procedure  must  be  performed  at  each  sta- 
tionary point,  and  then  the  results  are  added  to  obtain  the  final  answer. 

The  stationary  points  of  kR  occur  at  0 = 0 and  Tt , that  is,  where 


d ( k R > 
dO 


ka  sin  8 sin  8 sin  o 
o v 


= 0 . 


The  integral  over  the  expansion  of  kR  about  0 = 0 is 


2n 

J exP 

o 


ik  H + 

o 


iko2  R” 


2tt 


d0=  exp|ikRo]  J exp 


2 2 
ik  a sin  8 sin  8 0 

o Y 

2~~R 


1 71  sin  1 

(V) 

y 2 1 

f 2 ak  sin  6 sin  6 

o 

d0 


(3) 


where  R()  = 2a  sinZ  f(8  a~8  )/2]  and  the  asymptotic  limits  of  the  Fresnel  sine  and 
cosine  functions  were  utilized,  ^ 

X flradsliteyn,  I.S.,  and  Ryzkik,  1.  M.  (1965)  Table  of  Integrals  Series  and 
Products,  Academic  Press,  New  York  and  London,  pp  930-932. 
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Similarly  for  0 = v.  R^  = 2a  sin  |(0q  + 9)12]  and 


2tj 

J <-*xP 


ik©2  H ’ 


ik  H 


2jt 


d0  = expjikK^]  J exp 


2 2 
-ik a sin  0 sin  0 0 

o 


TIT 


a > > 1 


exp[ik  R^  - in 


'■>]  Vrf  sin  0 sin  0 ‘ 

* o 


d0 


(4) 


The  sum  of  the  contributions  from  the  two  stationary  points  yields 

> jfl  * l 

Vp,--Sx — 1 


where  1 - F-B;  here, 


(5) 


_ lx sin  0 

*■  ^07  exp 


i2ak  -in(^-)+-r|y-ta(^)  ('  ^ 1 ft 

\ 2 ak  sin  ^ ° 


arises  from  the  forward  moments  near  0=0  while 

I i 2 ak  sin 


B =Jl-  — 
j sin 


T oxp 


/vm  i.  rrvM /,  i 

y sin\  2 / [1  + pr-Tin 

V 2aksin(-^_j 


arises  from  the  backward  moments  near  0 = Jl.  Curves  for  |b|,  |f|  and  |l|  vs  0 
are  plotted  in  Figure  3. 


Figure  2.  Schematic  Representation 
of  the  Backward  (b)  and  Forward  (f) 
Fresnel  Zones  for  the  Magnetic  Field 
at  P Due  to  Current  Dipole  Moments 
Along  Ring  QQ'Q 
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Figure  3.  Plot  of  | I*'| , I B | , | b|  + | F 
| Bj  - | F | , and  | B - F | for  0 s 0 50 
Where  9 = 0.  1 radian 


Notice  that  as  long  as  the  near  field  terms  are  negligible,  the  contribution  from 


the  backward  moments  is  larger  than  from  the  forward  moments,  in  the  ratio 


mi 


(6) 


This  perhaps  surprising  result  arises  from  the  fact  that  the  backward  Fresnel 
zone  is  larger  than  the  forward  zone  (as  shown  in  Figure  21  in  a way  that  more  than 
compensates  for  the  additional  distance. 

The  apparent  contradiction  with  Huygen's  principle  is  resolved  in  a subsequent 

4 

9 - integration  over  all  such  rings  for  0 s 0 s 9 . This  integration  leads  to  a 
second  application  of  Kelvin's  principle  that  eliminates  the  backward  Fresnel  zone 
contribution. 

Ti  born,  M.  , and  Wolf,  E.  ( 1 9 7 01  Principles  of  Optics,  Pergamon  Press, 

Oxford  Fourth  Edition,  pp  370-375. 
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3.  NUMERICAL  ANALYSIS 


Eq.  (2)  can  be  numerically  integrated  via  t'hebyshev’s  quadratures,  which  is 
potentially  more  exact  than  the  Kelvin  approximation  because  the  moment  ring  can 
be  sampled  at  many  points.1’  Rewrite  Eq.  (2)  as 

i A 1 

VP)  = _ !c  ' (7) 


N exp  (i'j2  aka)  . , x 

a --  ojrtlT-  c 

n=  1 a ' \ n ' 


( 1 - cos  0 cos  0 


sin  0 sin  0 cos  0 ) 
o n 


sin  0 (cos  0 n - sin  0 sin  0 o - cos  0 cos  0 q cos  0n> 


0 = 7i  ( 2n  - 1 ) / 2 N . 
vn 

1 was  evaluated  on  a CDC  C6 00  computer  for  the  following  parameters: 
c 


= 30,000 
= 636  7.951km 


= 0.  1 rad. 


This  value  of  0 locates  the  observation  point  1’  at  a geodetic  distance  of 
636.8  km  from  pole  T in  f igure  1.  Table  1 gives  the  magnitudes  of  1 and  Ic>  and 
Table  2 gives  the  corresponding  phases.  For  0.  01  s fl  s 0.  09  there  appears  a ran- 
dom difference  in  the  magnitudes  of  1 and  1(„  that  is  less  than  0.  6 percent,  and  a 
random  difference  in  phase  that  is  less  than  0.  1 percent  (of  a maximum  possible  2jt ). 
This  phase  difference  corresponds  to  10  nsec  for  a 100  kHz  KORAN  signal. 

On  the  other  hand,  for  0 s 10~J  rad  (that  is,  for  rings  less  than  about  2A  it. 
radius)  1 and  I diverge.  An  approximation  of  the  exact  integral,  valid  for  small  0 , 
may  be  derived  from  Eq.  (2):  namely. 


5.  Abramowitz,  .VI.,  and  Stegun,  I.  A.  (1964)  Handbook  of  Mathematical  !■  unctions. 
National  Bureau  of  Standards  Applied  Maths.  Series,  55,  p 88P. 
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where  , 

I'  = tt  sin  0 exp  (i?))  J (7)  6 / 2)  [(2/tj)  + i (1  - 62)]  , 

H = a k y 2(  1 - cos  0 cos  0 Q)  , 

6 = sin  0 sin  0 /(I  - cos  0 cos  0 ) , 

o o 

and  Jj  (x)  is  Bessel's  function  of  first  order. 

This  expression  fo>-  i'  is  not  derivable  from  I in  Eq.  (5).  In  the  second  columns 
of  Tables  1 and  2,  the  values  of  l'  (marked  with  asterisks)  were  listed  for  0 s 10~4  rad, 
where  there  was  excellent  agreement  between  1 and  1 ' _ 


Table  1.  Comparison  of  Kelvin's  Approximation  |l|  With  Chebyshev's 
Quadratures  |lcl 


0 

\i\ 

III 

1 

o 

o 

c 

c c 

. 000001 

2. 09231  x 10'8V 

2. 09231  x 10'8 

0. 

1. 51515  x 10-5 

+70000. 

. 00001 

2.  08772  x 10‘6* 

2.  08772  x 10'6 

0. 

4.  16843  x 10'5 

+ 2000. 

. 0001 

1.  66259  x 10~4* 

1.  66136  x 10-4 

+ 0.  07 

1.  12897  x 10~4 

-32. 

. 001 

2.  18594  x 10-5 

4.  12690  x 10*6 

+430. 

. 01 

9.  87999  x 10'4 

9.  93645  x 10'4 

-0.  6 

. 02 

3.  02082  x 10'3 

3.  02006  x 10'3 

+ 0.  03 

. 03 

8.  50989  x 10~4 

8.  47396  x 10'3 

+ 0.4 

. 04 

3.73989  x 10'3 

3.74241  x 10'3 

-0.07 

. 05 

3.  63608  x 10'3 

3.  63226  x 10'3 

+ 0.  1 

. 06 

2.  34100  x 10'3 

2.  34601  x 10'3 

-0.  2 

.07 

5.  31362  x 10'3 

5.  31393  x 10'3 

-0.  006 

. 08 

3.42070  x 10'3 

3.  41286  x 10'3 

+ 0.  2 

.09 

4.40749  x 10*3 

4.  42000  x 10"3 

-0.3 

.099 

4. 98026  x 10'3 

5.  01428  x 10‘3 

-0.7 

. 0999 

4.71620  x 10 "3 

4.  82516  x 10"3 

-2.3 

. 09999 

4,  93795  x 10*3 

4.  82721  x 10‘3 

+ 2.3 

. 099999 

5.  25122  x 10‘3 

4.  59873  x 10-3 

+ 14. 

| l'|  was  used  for  these  values. 


Table  2.  Comparison  of  Kelvin's  Approximation  With 
Chebyshev's  Quadratures  L\c 


1 

0 

U 

Li 

C 

100(/_I  - /_lc)/2ff 

. 000001 

5.  82201* 

2.  68117 

5.  82125 

*.01 

-50. 

. 00001 

5.  82201* 

2.  68123 

5.  82125 

+.01 

-50. 

. 0001 

5.  82267* 

5. 82193 

5.  82143 

+ . 02 
+ . 008 

. 001 

2.  83853 

3. 66014 

-13. 

.01 

5.70793 

5.70768 

+.  004 

. 02 

5.77161 

5.77079 

+ .01 

. 03 

0.  110873 

0.  113704 

-.  05 

. 04 

2.  29917 

2. 29834 

+ .01 

. 05 

2.49314 

2.49241 

+.01 

. 06 

4.48320 

4.48406 

-.01 

. 07 

4.99256 

4.  99130 

+.  02 

. 08 

5.  56764 

5.  56849 

-.01 

. 09 

0.  685514 

0.  685329 

+ . 003 

. 099 

0.  600346 

0.  595509 

+ . 08 

.0999 

6.  20815 

6. 22394 

-.  3 

. 09999 

1.  15077 

1. 08156 

+ 1.  26 

. 099999 

1. 37070 

1.  18344 

+ 3.  00 

• [_1  was  used  for  these  values. 


4.  CONCLUDING  RKMAKKS 

From  the  above  it  may  be  concluded  that  Kelvin's  approximation  is  valid  for 
2A/a<’0  <Q  Q and  I for  0 <2A/a.  Near  0 = 0q  any  quantitative  comparison  is  over- 
whelmed by  the  divergence  caused  by  the  near-field  term. 
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